Journal of Atmospheric and Solar-Terrestrial Physics 115-116 (2014) 87-94 



ELSEVIER 


Contents lists available at ScienceDirect 

Journal of Atmospheric and Solar-Terrestrial Physics 

journal homepage: www.elsevier.com/locate/jastp 



Wavelength dependence of solar irradiance enhancement during 
X-class flares and its influence on the upper atmosphere 

Yanshi Huang a> * Arthur D. Richmond b , Yue Deng c , Phillip C. Chamberlin ', Liying Qian b , 
Stanley C. Solomon b , Raymond G. Roble b , Zuo Xiao e 

a Department of Electrical and Computer Engineering, University of New Mexico, NM, USA 
b High Altitude Observatory, National Center for Atmospheric Research, Boulder, CO, USA 
c Department of Physics, University of Texas at Arlington, TX, USA 
d Solar Physics Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD, USA 
e Department of Geophysics, Peking University, Beijing, China 



CrossMark 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 26 April 2013 
Received in revised form 
25 September 2013 
Accepted 11 October 2013 
Available online 5 November 2013 


Keywords: 
Solar flare 
F1SM 
TIE-GCM 


The wavelength dependence of solar irradiance enhancement during flare events is one of the important 
factors in determining how the Thermosphere-Ionosphere (T-I) system responds to flares. To investigate 
the wavelength dependence of flare enhancement, the Flare Irradiance Spectral Model (F1SM) was run for 
61 X-class flares. The absolute and the percentage increases of solar irradiance at flare peaks, compared 
to pre-flare conditions, have clear wavelength dependences. The 0-14 nm irradiance increases much 
more (~680% on average) than that in the 14-25 nm waveband ( — 65% on average), except at 24 nm 
(~220%). The average percentage increases for the 25-105 nm and 122-190 nm wavebands are ~120% 
and ~35%, respectively. The influence of 6 different wavebands (0-14 nm, 14-25 nm, 25-105 nm, 105- 
120 nm, 121.56 nm, and 122-175 nm) on the thermosphere was examined for the October 28th, 2003 
flare (X17-class) event by coupling FISM with the National Center for Atmospheric Research (NCAR) 
Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIE-GCM) under geomagneti- 
cally quiet conditions (I(p= 1). While the enhancement in the 0-14 nm waveband caused the largest 
enhancement of the globally integrated solar heating, the impact of solar irradiance enhancement on the 
thermosphere at 400 km is largest for the 25-105 nm waveband (EUV), which accounts for about 33 K of 
the total 45 K temperature enhancement, and ~7.4% of the total ~ 11.5% neutral density enhancement. 
The effect of 122-175 nm flare radiation on the thermosphere is rather small. The study also illustrates 
that the high-altitude thermospheric response to the flare radiation at 0-175 nm is almost a linear 
combination of the responses to the individual wavebands. The upper thermospheric temperature and 
density enhancements peaked 3-5 h after the maximum flare radiation. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

A solar flare results from a sudden, intense release of magnetic 
energy in the atmosphere of the Sun, which produces rapid 
increases in electromagnetic radiation, from gamma rays to radio 
wavelengths (Einar and Gordon, 1988). Flares are classified as A, B, 
C, M, or X according to the maximum flux of 0.1 -0.8 nm X-rays 
measured near the Earth (Garcia, 2000). While C-class flares are a 
common occurrence during years near solar maximum, the 
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frequency of X-class flares during solar maximum is approximately 
two per month (http://spacemath.gsfc.nasa.gov). There are large 
spectral variations among flares (Thomson et al., 2004; Tsurutani 
et al., 2005; Chamberlin et al., 2012). The enhancement of the 
Extreme Ultraviolet (EUV, ~25 to — 120 nm) spectral irradiance 
depends on the location of a flare, whereas the flare enhancement 
of X-ray Ultraviolet (XUV, ~0.1 to ~25nm) radiation depends 
only weakly on location (Qian et al., 2010). 

Solar EUV photons are the primary energy source of the neutral 
and ionized constituents of the Thermosphere-Ionosphere (T-l) 
system (Mitra, 1974; Liu et al., 2011). The extra ionization in the 
ionosphere caused by flares increases electron density, which 
influences the absorption and refraction of radio waves propagat- 
ing through the ionosphere from one station to another (Einar and 
Gordon, 1988). Also, flares are often associated with Coronal Mass 
Ejections (CMEs), which may cause significant geomagnetic storms 
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(Munro et al., 1979; Harrison, 1987; Tsurutani et al., 1988; 
Gonzalez et al., 1994; Michalek, 2009; Bein et al., 2012). Previous 
studies have shown that the impact of flares on the T-I system 
varies because flares have different magnitudes, locations on the 
solar disk, rise rates and decay rates (Tsurutani et al., 2005; Sutton 
et al., 2006; Qian et al., 2010, 2011; Zhang et al., 2011, 2012). Since 
XUV dominates ionization in the lower thermosphere ( < 150 km), 
and EUV dominates in the upper thermosphere (Qian et al., 2011 ), 
the impact of solar flares on the upper atmosphere depends on the 
solar irradiance enhancement at different wavelengths. 

The purpose of this paper is to investigate how different wave- 
bands of solar irradiance are enhanced during flares and how they 
impact the global thermosphere. The Flare Irradiance Spectral Model 
(FISM) (Chamberlin et al., 2007, 2008) was employed to estimate the 
spectra of solar irradiance. The National Center for Atmospheric 
Research (NCAR) Thermosphere-Ionosphere-Electrodynamics General 
Circulation Model (TIE-GCM) (Roble et al., 1988; Richmond et al., 1992) 
driven by the FISM was used to simulate the thermospheric responses 
to flares. 

2. Model description 

2.1. Flare irradiance spectral model (FISM) 

FISM is an empirical model of the solar irradiance spectrum 
from 0.1 to 190 nm at 1 nm resolution and on 1 -minute time 
cadence. The high temporal resolution of FISM makes it possible to 
study the variations due to solar flares. This model is based on the 
data provided by the Solar Extreme ultraviolet Experiment (SEE) 
(Woods et al., 2005) on the Thermosphere Ionosphere Mesosphere 
Energetics and Dynamics (TIMED) satellite, and the SOLar STellar 
Irradiance Comparison Experiment (SOLSTICE) (Rottman, 2000) on 
the Solar Radiation and Climate Experiment (SORCE). 

FISM estimates the daily component of irradiance, including 
the variations from the solar cycle and solar rotation of active 
regions (Chamberlin et al., 2007). The flare component of FISM 
includes both the impulsive and gradual phase variations, and is 
based on a reference set of 39 large flares from 2002 to 2005 
measured by SEE and SORCE (Chamberlin et al., 2008). This 
reference set of measurements is fit to a flare proxy to determine 
the coefficients. The Geostationary Operational Environmental 
Satellites (GOES) 0.1 -0.8 nm fluxes are used as the flare proxy to 
empirically model the flare variation for FISM, because of its high 
temporal resolution and reliable data since 1970 with few data 
gaps, as well as plans for continued future measurements. A linear 
relation is found between the irradiance of GOES 0.1 -0.8 nm and 
the irradiance for the X-rays at wavelengths less than 14 nm, while 
a power law relation with exponent of 0.647 is found between the 
irradiance of GOES 0.1 -0.8 nm and the irradiance for EUV at 
wavelengths larger than 14 nm (Chamberlin et al., 2008). 

Due to large variations during flares for different wavelengths, 
and due to limited measurements, the FISM flare uncertainty has a 
wavelength dependence that varies from 10% to above 100%. 
However, compared to models with only daily outputs, FISM 
improves the estimation of solar flares significantly. More flare 
data can help reduce the FISM flare uncertainties. The newly 
available and more accurate data from the Extreme ultraviolet 
Variability Experiment (EVE) on the Solar Dynamics Observatory 
(SDO) (Woods et al., 2010, 2011) will add the additional variable of 
temporal delays in some of the EUV emissions during the gradual 
phase of flares (Chamberlin et al., 2012). 

2.2. NCAR TIE-GCM 

The NCAR TIE-GCM is a first-principles, three-dimensional, non- 
linear representation of the coupled thermosphere and ionosphere 


system. It solves the momentum, energy and continuity equations for 
neutral and ion species in pressure coordinates (Roble et al., 1988), 
with a self-consistent calculation of ionospheric wind dynamo effects 
(Richmond et al., 1992). The external forcings of TIE-GCM are mainly 
the solar irradiance, magnetospheric energy, and tidal perturbations 
at the lower boundary of the model. Magnetospheric energy inputs 
include auroral particle precipitation and high-latitude ion convec- 
tion. In this study, the Heelis potential model (Heelis et al., 1982) was 
used to specify the high-latitude electric field while the auroral 
particle precipitation is as described by Roble and Ridley (1987). The 
latest version (vl.94) of the TIE-GCM was run with a 5° x 5 x half 
scale height resolution (longitude x latitude x vertical). 

Solar irradiance spectra estimated by the FISM were used as 
solar input to the TIE-GCM in our study. The spectral range 
required by the TIE-GCM is 0.05-175 nm, so the part of FISM 
spectra at wavelengths greater than 175 nm is not included. The 
low-resolution binning scheme developed by Solomon and Qian 
(2005) is utilized to separate the FISM spectrum into 37 bins for 
the solar EUV energy deposition calculation. The spectrum 
between 0.05-105 nm is divided into 22 non-uniform bins accord- 
ing to the changes of cross sections and photon energy within a 
bin. The part of spectra with wavelengths longer than 105 nm is 
divided evenly into 5-nm bins except for the Lyman— a line 
(121.56 nm) (Qian, 2007). 

The heating due to solar irradiance from 0.05 to 105 nm is 
calculated through integrating heating over the first 22 bins as 
shown by Wang (1998) in Eq. (2.8), which depends on the solar 
flux at different wavelengths, the ionization cross sections of 
N 2 ,0 2 , and O, and the column number densities of N 2 ,0 2 , and O. 
The contribution of the Schumann-Runge continuum (105- 
175 nm) is integrated over the remaining 15 bins. The photodisso- 
ciation by radiation in the Schumann-Runge continuum is the 
primary direct dissociation process for 0 2 in the thermosphere. 
The excess energy from 0 2 dissociation goes directly into heating, 
which is about 6.6% of the absorbed Schumann-Runge radiation 
(Schunk and Nagy, 2009). The quenching process of 0( 1 D) resulted 
in the 0 2 dissociation also contributes to the neutral heating. The 
neutral gas heating by the Schumann-Runge bands is specified by 
the empirical formula in Strobel (1978), with a heating efficiency 
of about 0.3. 


3. Results and discussion 

3.1. Solar irradiance enhancements from the FISM for X-class flares 

A flare is classified as A, B, C, M or X according to the peak flux 
of 0.1 -0.8 nm irradiance near the Earth. Within each class, a linear 
scale from 1.0 to 9.9 is used to specify the magnitude of a flare, 
except for X-class flares, the class number can exceed 9.9. The 
biggest flare in the National Oceanic and Atmospheric Adminis- 
tration (NOAA) records is an X28 flare that occurred on November 

4, 2003. The November 4, 2003, flare saturated GOES X-ray sensor 
(XRS), and its magnitude was initially estimated as X28 through 
extrapolating the GOES curve (http://www.swpc.noaa.gov). 
Thomson et al. (2004) argued that the magnitude of this X28 flare 
should be as large as X45 5 based on the measurements of 
ionospheric VLF (Very Low Frequency) radio phase change. 

Fig. 1 depicts the FISM outputs for the X17 flare on October 
28th, 2003 (day of year 301). The top panel shows the temporal 
variation of total solar flux integrated over 0-190 nm wavelengths, 
which started to increase at 11 :00 UT and reached a maximum flux 
of 100 mW/m 2 at around 11:07 UT. Solar spectra before the flare 
and at the flare peak are shown in the middle panel. The black line 
is for the pre-flare time which is marked with a triangle in the top 
panel, whereas the red line is for the time at the flare peak which 
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Fig. 1. F1SM solar irradiance results for the X17 flare on October 28th, 2003 (day of 
year 301). Top: time variation of total solar irradiance integrated over the 
wavelengths from 0 to 190 nm. The flare reached a maximum flux of 0.1 W/m 2 
at around 11:07 UT. Middle: solar spectra before the flare at the moment denoted 
by the triangle in the top panel (black) and at the flare peak denoted by the 
diamond (red). Bottom: percentage increase of solar irradiance comparing the peak 
and pre-flare conditions. (For interpretation of the references to color in this figure 
caption, the reader is referred to the web version of this article.) 


is marked with a diamond in the top panel. The clear peak at about 
121.56 nm before the flare is the emission line of hydrogen Lyman- 
a, which will be discussed in more details later. The flux enhance- 
ment varied with wavelength, as illustrated in the bottom panel, 
which shows the percentage increase of solar flux at the flare peak 
compared to the pre-flare condition. The solar irradiance increased 
the most in the 0-14 nm range, about 1500% on average, whereas 
it only increased ~100% in the 14-25 nm range, except at 24 nm 
(~300%). The average enhancements of irradiance in the 25- 
105 nm and 122-190 nm wavebands are about 230% and 70% for 
this flare, respectively. The solar irradiance increase modeled by 
the FISM is for the brief, initial enhancement during the impulsive 
phase of this flare. 

In order to have a better understanding of the wavelength 
dependence of irradiance enhancement of X-class flares, we 
examined the spectra of solar irradiance percentage increase 
during 61 X-class flares between 1989 and 2012 using the FISM. 
These flares ranged from X1.0 to X28, and mostly occurred in the 
maximum and descending phases of solar cycles. The 20 X-class 
flares occurred during 2002-2005 among these 61 flare events 
along with another 19 M-class flares are the 39 large flares that are 
used as the basis of FISM (Chamberlin et al., 2008). Fig. 2 shows 
the percentage increase of the 61 X-class flares. The red line 
denotes the October 28th, 2003 flare event, and the green line 
is the average of the 61 X-class flares. The 0-14 nm irradi- 
ance increases much more (~680% on average) than that for the 
14-25 nm waveband ( ~ 65% on average), except for the peak at 
around 24 nm (~220%). The average percentage increases for 
the 25-105 nm and 122-190 nm wavebands are ~120% and 
~ 35%, respectively. There is a relatively small enhancement in 
the 14-25 nm wavelength range, which contains many strong 
solar emission lines, such as Fe IX. When a flare occurs, these 
emissions do not increase very much, and sometimes even 
decrease (Woods et al., 2011). This is because most of these 
emissions are formed in the solar corona at the approximate 
temperature of the corona, around 1-2 million K. When the source 
ions are heated to much higher temperatures, Fe IX quickly ionizes 
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Fig. 2. FISM solar irradiance results for the percentage increases comparing the 
peak and pre-flare conditions for the 61 X-class flares that occurred between year 
1989 and 2012. The red line denotes the October 28th, 2003, flare event, and the 
green line is the average of the 61 X-class flares. The low enhancement at 15-25 nm 
is a result of the depletion of source ions, whereas the low enhancement at 115— 
129 nm, except for the 121.56 nm bin, is rather instrumental. The lower envelope 
longer than 25 nm is for the X1.5 flare that occurred on April 21, 2002. (For 
interpretation of the references to color in this figure caption, the reader is referred 
to the web version of this article.) 


to higher states, depleting the source ion population for Fe IX 
emissions, so that there is a relatively small enhancement in the 
15-25 nm wavelength range. On the other hand, the low enhance- 
ment in the 115-129 nm range, excluding the 1 nm hydrogen 
Lyman-a bin at 121.56 nm, is rather instrumental. In the TIMED/ 
SEE instrument, which FISM is based on, an aluminum filter had to 
be put in place to block 99% of the Lyman-a emission line at 
121.56 nm so that it would not saturate the detector and still be 
enough to be measured (Eparvier et al., 2001 ). This filter blocks the 
entire 115-129 nm band, which is then assumed to vary in the 
same manner as the Lyman-a emission that is measured by SEE. 

In this study, the XUV/EUV spectrum is divided into six 
wavebands according to the penetration depth in the atmosphere 
and the irradiance enhancement during flares: 0-14 nm and 
14-25 nm wavebands for XUV, 25-105 nm and 105-120 nm for 
EUV, the Lyman-a line (121.56 nm), and the 122-175 nm wave- 
band for far ultraviolet (FUV) irradiance. We will investigate how 
these wavebands contribute to thermosphere responses to the 
X17-class flare that occurred on October 28th, 2003, one of the 
most severe flares in the last solar cycle. 

3.2. Temporal variation of thermosphere response at 400 km 

To evaluate the influence on the thermosphere of the 0-14 nm 
radiation for the October 28th, 2003 flare, the TIE-GCM was run 
holding all radiation constant at the pre-flare level except for that 
in the 0-14 nm band, which was varied as prescribed by the FISM. 
The difference between this run and a base TIE-GCM run with 
constant radiation in all wavebands was then evaluated. The 
influence of solar irradiance enhancement for the other wave- 
bands was calculated in the same way. To minimize geomagnetic 
influences, all the simulations are under geomagnetically quiet 
conditions (Kp= 1). In this study, the globally integrated or 
averaged values were analyzed to focus on large-scale responses. 

The temporal variations of solar irradiance flux in the 0-14 nm, 
25-105 nm and 122-175 nm wavebands were selected and empha- 
sized, since they contain most of the solar irradiance enhancement 
(Fig. 3). The solar flux variations for different wavebands calculated by 
the FISM are depicted in the first row. The solar flux in the 0-14 nm 
waveband is quite small before the flare ( ~ 1 mW/m 2 ), as compared 
with that in the 122-175 nm waveband ( ~ 12 mW/m 2 ), but it 
increased the most ( ~ 20 mW/m 2 ) at the flare peak. It is also noticed 
that the peak time of solar flux varies. Although the total flux 
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Fig. 3. Time variations of different wavebands for DOY 301, first column for 0-14 nm, second column for 25-105 nm and third column for 122-175 nm. Top row: solar 
irradiance calculated by the FISM. Second row: perturbation of global heating due to solar irradiance. The time integration (one day) of global solar heating (TIGSF1) deposited 
is labeled in red. Third row: perturbation of globally integrated Joule heating. Fourth row: Globally averaged perturbation of temperature at 400 km. Bottom row: Globally 
averaged perturbation of neutral density at 400 km. (For interpretation of the references to color in this figure caption, the reader is referred to the web version of this 
article.) 


integrated from 0 to 190 nm peaked at 11 :07 UT, the fluxes for 
105-120 nm and 122-190 nm peaked earlier at 11:05 UT, those for 
25-105 nm and 121 nm peaked at 11:06 UT, the flux for 14-25 nm 
peaked at 11:07 UT, and that for 0-14 nm peaked at 11 :10UT. 
Therefore, except for the Lyman-a emission line, basically the peak 
time decreases with the increasing wavelength. The second row 
shows the enhancement of globally integrated solar heating in the 
upper atmosphere during the flare, which was calculated by taking 
the difference between the flare and non-flare runs for each wave- 
band. The peak solar heating for the 0-14 nm waveband is the largest 
(~ 1216 GW), which is more than three times as large as the peak for 
the 25-105 nm waveband (~ 379 GW). The time integration (one 
day) of global solar heating (TIGSH) for the 25-105 nm waveband 
(1.85xl0 10 J) is only half of that for the 0-14 nm waveband 
(~3.48 x 10 10 J) but is larger than that for the 122-175 nm wave- 
band ( ~ 1 .08 x 10 10 J). In the third row, Joule heating enhancement 
due to increased Pedersen conductivity shows quite different 
responses in time and magnitude for different wavebands. There is 
a rather rapid and large enhancement of Joule heating responding to 
the solar flare irradiance for the 0-14 nm waveband. However, for 
122-175 nm, the peak of the Joule heating enhancement has an 
obvious delay with respect to the flare peak. The response for 25- 
105 nm is a combination of the patterns for the 0-14 nm and the 
122-175 nm wavebands. The radiation at 122-175 nm impacts the 
thermosphere through dissociation and heating, but not ionization, 


and is absorbed by 0 2 , which has significant densities only in the 
lower thermosphere. Therefore, it affects the Pedersen conductivity 
indirectly through affecting the neutral temperature and density in 
the lower thermosphere, which in turn affect the ion density and ion- 
neutral collision frequency. Due to the fact that the neutral tempera- 
ture and density changes persist longer than the ionization enhance- 
ments, the associated changes of Pedersen conductivity and Joule 
heating also persisted hours after the flare. The global Joule heating 
increased, but the magnitude was quite small compared with the 
global heating due to the solar irradiance. Joule heating in these 
results was rather small because the imposed high-latitude electric 
field was relatively small. 

The fourth and bottom rows of Fig. 3 depict the thermospheric 
temperature and density perturbations at 400 km, which is close 
to altitudes of some Low-Earth Orbiting (LEO) satellites, such as 
the Challenging Minisatellite Payload (CHAMP). The total tem- 
perature and neutral density perturbations at 400 km were ~45 K 
and — 11.5%, respectively. The thermospheric perturbations at 
400 km due to solar irradiance enhancement at 25-105 nm 
(~33 K increase of temperature: ~7.4% increase of neutral den- 
sity) are much larger than at 0-14 nm ( ~8 K increase of tempera- 
ture; ~2.5% increase of neutral density). The enhancement of the 
122-175 nm waveband resulted in little effect on the high-altitude 
thermospheric responses. Most of the irradiance in 122-175 nm is 
absorbed by 0 2 in the lower thermosphere. Some of the heat 
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Table 1 

Heating efficiency of solar irradiance at different wavelengths to thermosphere at 
400 km for the X17 flare event on October 28th, 2003 (day of year 301 ). 



0-14 nm 

25-105 nm 

122-175 nm 

S SH peak (GW) 

1216 

379 

548 

TIGSH (J) 

3.48el0 

1.85el0 

1.08el0 

<5 Tn peak (K) 

8.4 

32.8 

1.2 

(5/iPeak 

2.50% 

7.39% 

0.58% 

S Tn/S SH (K/GW) 

0.0069 

0.0865 

0.0022 

S Tn/TIGSH (K/J) 

2.41 e- 10 

17.73e- 10 

1.09e — 10 

ap/«SH (%/GW) 

0.0021 

0.0200 

0.0011 

S^/TIGSH (%/J) 

0.7e — 10 

4.0e- 10 

0.5e- 10 


deposited in the lower thermosphere conducts up to the upper 
thermosphere and produces small but long-lasting thermospheric 
perturbations. On the other hand, the 0-14 nm irradiance is 
absorbed by O and N 2 as well as 0 2 . Therefore , it has a larger 
fraction absorbed in the upper thermosphere compared to the 
Schumann-Runge continuum (Torr et al., 1979; Roble and Emery, 
1983). As the high-altitude atmosphere conducts heat downward 
over time, the upper thermosphere cools down. 

Although the largest heating enhancement comes from solar 
irradiance in the 0-14 nm wavelength range, most of the thermo- 
spheric perturbations at 400 km are due to the 25-105 nm 
waveband (Table 1). The heating efficiency, which is defined as 
temperature (density) perturbations divided by the peak solar 
heating, is an order larger for the 25-105 nm waveband than those 
for the other two wavebands. Therefore, the solar irradiance in the 
25-105 nm wavelength range influences the thermosphere at 
400 km most effectively. The maximum high-altitude thermo- 
spheric responses occurred 3 ~ 5 h later than the flare peak. 

Fig. 4 compares the full simulation including all the variations 
in the 0-175 nm wavelength range with the summation of 
perturbations due to the six individual wavebands (0-14 nm, 
14-25 nm, 25-105 nm, 105-120 nm, 121.56 nm, and 122-175 nm). 
The black lines represent the globally averaged perturbations of 
(a) solar heating deposited into the thermosphere, (b) neutral 
temperature and (c) density at 400 km due to the whole flare 
spectrum from 0 to 175 nm, while the red lines represent the 
summations of perturbations at 400 km due to the six separate 
bands. The black lines and red lines are very close to each other, 
which illustrates the nearly linear response of the thermosphere to 
different wavebands of solar irradiance at 400 km. Therefore, the 
high-altitude thermospheric response to the 0-175 nm flare radia- 
tion is almost a linear combination of the responses to individual 
wavebands. This suggests a possibility to predict the variations of 
temperature and neutral density at LEO due to flares using the 
solar irradiance measurements at different wavebands. For future 
flares, measurements covering the different wavebands of solar 
irradiance can give estimations of the thermospheric responses 
separately, and then the total response is predicted as a linear 
combination of them. 

3.3. Altitudinal distribution of the thermospheric disturbance 

The results for the 0-14 nm and 25-105 nm wavebands are 
further discussed in this section since they have the most 
significant influence on the thermosphere. Fig. 5(a) depicts the 
altitudinal distributions of heat perturbations due to solar irra- 
diance enhancement in different wavebands at the flare peak. The 
largest perturbations of heating per unit mass occur above 150 km, 
which is dominated by the 25-105 nm waveband, while the 
heating below 150 km is dominated by the irradiance of 0-14 nm. 
The solar heating rate per unit mass peaks in the upper thermo- 
sphere around 300 km at high solar activity conditions, and it 


a 



Time 


b 



Time 


c 



Time 


Fig. 4. Comparison between the full simulation that includes the flare enhance- 
ments of all wavebands, and the summation of perturbations due to the 6 individual 
bands (0-14 nm, 14-25 nm, 25-105 nm, 105-120 nm, 121.56 nm, and 122-175 nm). 
The black lines show the variations using the full 0-175 nm FISM solar spectra, 
while the red lines represent the sum of variations due to the FISM spectra in the 
separate wavebands, (a) Globally averaged solar heating perturbation at 400 km. 
(b) Globally averaged temperature perturbation at 400 km. (c) Globally averaged 
neutral density perturbation at 400 km. (For interpretation of the references to 
color in this figure caption, the reader is referred to the web version of this article.) 

is dominated by the heating process of neutral-neutral and ion- 
neutral chemical reactions (Roble et al., 1987). The heating 
perturbation resulted in similar temperature response, as shown 
in Fig. 5(b). The temperature perturbation below ~ 150 km is 
larger for the solar heating in the 0-14 nm waveband, while for 
the upper thermosphere above 150 km, the perturbation is domi- 
nated by the heating at 25-105 nm. Consequently, the neutral 
density perturbation around 400 km is dominated by the solar 
heating in the 25-105 nm waveband as illustrated in Fig. 5(c). 
Therefore, the solar irradiance at 25-105 nm is significant for high 
altitude heating at the flare peak, while the irradiance at 0-14 nm 
plays a dominant role below 150 km. 

The temporal and altitudinal responses of the thermosphere to 
solar flare irradiance varied significantly with wavelength. Fig. 6 
shows the temporal variations of the altitudinal distributions for 
different wavebands: the first column for the 0-14 nm waveband, 
the second column for 25-105 nm and the third column for the full 
waveband from 0 to 175 nm. As shown in the first row, the largest 
percentage increase of heating due to solar irradiance in the 0-14 nm 
waveband is below 150 km, while that for 25-105 nm is above 
150 km. The total heating due to solar irradiance for the full waveband 
(0-175 nm) increased up to 140% below ~ 150 km. The second and 
third rows show the corresponding perturbations of temperature and 
neutral density, respectively. The maximum temperature perturbation 
for 0-175 nm was about 47 I< at high altitudes, with ~ 9 I< contributed 
by the 0-14 nm waveband and ~34 K by the 25-105 nm waveband. 


92 


V. Huang et al. / Journal of Atmospheric and Solar-Terrestrial Physics 115-116 (2014) 87-94 


a 



b 



c 



Fig. 5. Altitudinal distribution of TIE-GCM simulation results at the flare peak due 
to solar irradiance enhancement in different wavelength ranges, (a) Globally 
averaged heating perturbation, (b) Globally averaged temperature perturbation, 
(c) Globally averaged neutral density perturbation. The density perturbations 
below 145 km can sometimes be negative, and only positive values are plotted. 


Similarly, the maximum neutral density enhancement due to the 
0-175 nm waveband was around 18.8%, with ~3.7% contribution 
from 0-14 nm and ~12.9% from 25-105 nm. Apparently, the solar 
irradiance at 25-105 nm plays a dominant role in impacting the upper 
atmosphere. It can be explained by the fact that the upper thermo- 
spheric response is much stronger to the heating in the upper 
thermosphere than to the heating in the lower thermosphere (Deng 
et al., 2011; Huang et al., 2012). 

The white lines in the middle and bottom rows of Fig. 6 
represent the moments of the largest temperature and density 
enhancements at different altitudes. The moments of maximum 
temperature perturbations denoted by the white line (second row) 
correspond well to the solar heating at different altitudes (first 
row). For the 0-14 nm waveband, the maximum temperature 
enhancement at about 150 km occurs quickly corresponding to 
the maximum heating at the same altitude. Then the heat 
conducts upward and downward, gradually heating the surround- 
ing thermosphere. For the 25-105 nm waveband, since most of the 
heat is deposited above 150 km, the heat conduction is mainly 
downward. Therefore it takes longer time to heat the lower 
thermosphere. The peak upper thermospheric perturbation had 
3-5 h time delay compared to the flare peak. The large-scale 
gravity waves, which are generated by thermospheric heating 
and expansion on the dayside, transport energy from dayside to 
nightside. Pawlowski and Ridley (2008) investigated the global 
thermospheric response to the October 28, 2003, flare and found 
that the gravity wave propagates at the speed of local sound speed 
plus the neutral wind velocity. The global mean sound speed at 
400 km is around 900 m/s and the average neutral wind speed is 
around 200 m/s. Therefore, with a total speed of 1100 m/s, it takes 
about 5 h for the wave to travel from dayside to nightside. 

For this study, in order to focus on the solar flare influence on the 
upper atmosphere and minimize the geomagnetic effect, we did a 
non-realistic simulation with constant I<p= 1. Therefore, the 


0-14 nm 


25-105 nm 


0- 1 75 nm 





Fig. 6. Temporal variations of the altitudinal distribution of TIE-GCM simulation results due to solar irradiance enhancement in the 0-14 nm, 25-105 nm and 0-175 nm 
wavelength ranges. First row: globally averaged percentage increase of solar heating that deposited in the thermosphere. Second row: globally averaged temperature 
enhancement. Third row: globally averaged percentage increase of neutral density. The white lines represent the moments of the largest temperature and density 
enhancements at different altitudes. The density enhancement below 145 km can sometimes be negative. 
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comparison between neutral density from our TIE-GCM runs and the 
observations from CHAllenging Minisatellite Payload (CHAMP)/GRACE 
Gravity Recovery And Climate Experiment (GRACE) could not be used 
to validate our model results. 


4. Conclusions 

The responses of the thermosphere to solar flares vary in 
magnitude and duration, since flares have different energy deposi- 
tion heights, magnitudes and durations for different wavebands. 
The wavelength dependence of the solar flare enhancement is one 
of the important factors in determining how the terrestrial atmo- 
sphere responds to flares. The FISM spectra show that for the 61 X- 
class solar flares that occurred between 1989 and 2012, the solar 
irradiance enhancement at the flare peak relative to the pre-flare 
condition has a clear wavelength dependence. The 0-14 nm 
irradiance increases much more (~680% on average) than that 
for the 14-25 nm waveband (~65% on average), except at around 
24 nm (~220%). The average percentage increases for the 25-105 nm 
and 122-190 nm wavebands are ~120% and ~35%, respectively. 

The thermospheric response to six different wavebands of solar 
irradiance for the X17-class solar flare on October 28th, 2003, was 
investigated through coupling the FISM to the TIE-GCM as solar 
input. The globally integrated solar heating and the one-day 
integration of it are largest in the 0-14 nm waveband. However, 
the impact of solar irradiance enhancement on the upper thermo- 
sphere is largest for the 25-105 nm waveband (EUV). At ~400 km, 
it accounts for about 33 K of the total 45 K temperature enhance- 
ment, and ~7.4% of the total — 11.5% neutral density enhance- 
ment. For this particular simulation, the maximum temperature 
perturbation is ~47I< (34 K due to EUV irradiance), and the 
maximum neutral density perturbation is ~18.8% (12.9% due to 
EUV irradiance) at higher altitude than 400 km. The effect of 122- 
175 nm (FUV) irradiance is small. The high-altitude temperature 
and neutral density enhancements maximize 3 - 5 h after the flare 
peak. The comparison between a full simulation including all the 
variation at 0-175 nm wavelengths and the summation of pertur- 
bations due to the six individual wavebands indicates that at 
400 km, there is a nearly linear response of the thermosphere to 
different wavebands of solar irradiance. It suggests a possibility to 
predict the variations of temperature and neutral density at LEO 
due to flares using the solar irradiance measurements at different 
wavebands. 
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